Solving the two-center nuclear shell-model problem with arbitrarily-orientated 

deformed potentials 



Alexis Diaz- Torres 

Department of Nuclear Physics, Research School of Physical Sciences and Engineering, 
Australian National University, Canberra, ACT 0200, Australia 
(Dated: October 2, 2008) 

A general new technique to solve the two-center problem with arbitrarily-orientated deformed 
realistic potentials is demonstrated, which is based on the powerful potential separable expansion 
method. As an example, molecular single-particle spectra for 12 C + 12 C — > 24 Mg are calculated using 
deformed Woods-Saxon potentials. These clearly show that non-axial symmetric configurations play 
a crucial role in molecular resonances observed in reaction processes for this system at low energy. 
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The description of a particle moving in the field of two 
fixed potential centers separated by a distance R (the 
two-center problem) is fundamental in classical [l[ and 
quantum mechanics 0, Q , and finds myriad applications 
in celestial mechanics, quantum chemistry, atomic and 
molecular physics and nuclear physics. This appears (i) 
in the study of the scattering of radiation by two black 
holes [3] , (ii)_} n the quantum mechanical theory of chemi- 
cal binding [5| , (hi) in the description of electron-positron 
pair production in heavy ion and ion-atom collisions Q, 
(iv) in the study of the properties of baryons contain- 
ing two heavy quarks (QQq) H], and (v) in phenomena 
related to nuclear molecules [8| . 

The applications to phenomena in low energy nu- 
clear physics [1] were first introduced (in practice) by 
the Frankfurt school using the two-center shell model 
(TCSM) based on a double oscillator potential @. Im- 
proved versions of this approach have been suggested for 
dealing with super-asymmetric fission (Io| and asymmet- 
ric fission with deformed fragments ll| . In all these mod- 
els, the two-center potentials are rotationally symmetric 
about the internuclear axis. These potentials are appro- 
priate, e.g., for the description of binary fission where the 
fragments are spherical or deformed with their intrinsic 
symmetry axis aligned with the internuclear axis. 

All alignments are possible in collisions of deformed 
nuclei. The major role of orientation of the deformed 
target in the onset of quasi-fission, whose understanding 
is very important to un puz zle the formation mechanism 
of superheavy elements [12J, has been demonstrated [13| 
by fission measurements in reactions forming heavy ele- 
ments. A more general TCSM is required for a proper 
description of these reactions within the molecular pic- 
ture [8j. It is justified at low incident energies near the 
Coulomb barrier, as the radial motion of the nuclei is 
expected to be adiabatically slow compared to the rear- 
rangement of the two-center mean field of nucleons. To 
my knowledge, only one attempt has been made to ac- 
count for arbitrarily-orientated deformed nuclei, in which 
the wave function expansion method (usual diagonaliza- 
tion procedure) and specific potentials (two ellipsoidally 





















y r 2 ™ 







FIG. 1: Schematic picture of the coordinates used to define 
the two-center potential |T]) in the collision between two de- 
formed nuclei. See text for further details. 



deformed Gausssian potentials) were applied to describe 
the reaction 13 C + 16 O [H|. 

I present a general new technique to solve the two- 
center problem with arbitrarily-orientated deformed frag- 
ments. The procedure is based on the powerful potential 
separable expansion method [l5j that has been success- 
fully used to solve the two-center problem with spherical 
Woods-Saxon (WS) potentials HEE^l- The technique is 
shown using deformed WS potentials, but it can also be 
employed with other types of deformed potential, pro- 
vided a suitable set of basis functions is selected. Hence, 
the method has applications in many areas of science. 
The formalism is described first, and illustrated after- 
wards with calculations of molecular single-particle (sp) 
spectra for the reaction 12 C + 12 C, which is of great as- 
trophysical interest [l8| . These calculations show the ma- 
jor role of non-axial symmetric configurations in forming 
a nuclear molecule, as the overlapping nuclei keep their 
identity and can "dance" for quite some time at contact. 

The finite depth nuclear potential V s (y,(3\q) of each 
spheroidal nucleus (s = 1,2) is chosen to be a de- 
formed WS with a spin-orbit term and deformation 
parameters (3\p. For protons, the Coulomb potential 



^Coui( r '^o) |20| of a uniformly charged spheroid with 
charge Z s e (Z s being the total charge of each fragment) 
should be added to the nuclear potential. This treatment 
of the Coulomb interaction is suitable for separated frag- 
ments, but may not be the most appropriate prescription 
for overlapping nuclei. Its validity could be assessed by 
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comparing this two-center Coulomb field to the field gen- 
erated by a uniformly distributed charge within a dinu- 
clear shape. The Coulomb interaction shifts up the pro- 
ton levels with respect to the neutron levels, and increases 
the sp potential barrier between the fragments 171 ] . 

The origin of the total deformed potential V s (r,f3xo) 
is placed at the position R s in the overall center-of-mass 
(CM) system, where its intrinsic symmetry axis is orien- 
tated by Eulerian angles fl s = (<j) s ,9 s ,0) (0 < 4> s < 2tt 
and < 9 S < n) with respect to the initial internuclear 
axis (see Fig[T]). Thus, the two-center potential is 
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where k — h~ l p is the sp wave-number operator and U 
is the operator of finite rotations 21| , whose inclusion is 
the key new aspect of this work. Of course, the rotation 
operator is not required for spherical nuclei fl7| . Each 
potential V s in |T]) is expressed as 
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within a truncated spherical sp harmonic oscillator basis, 
v = 1, . . . N}, with the spin- angular part having the 
total angular momentum j with projection m, e.g., in the 
momentum representation (see for further details) 

\u) = \nljm) = g nl (k) ■ [r^(k) ® X i («)]£,. (3) 

The number N of basis states is defined by l ma x (num- 
ber of partial waves in which the potential acts) and n max 
(the number of separable terms in each partial wave). 
The values of l ma x and n max are determined by the con- 
vergence of the sp energies, which is accelerated using 
the technique of the Lanczos c-factors (2^]. These are 
nmax=3 and Z maa; =4 for the studied reaction 12 C + 12 C. 

For bound (or quasi-stationary) states, the formal so- 
lution of the Schrodinger equation is 



\<p) = G (E)V\<p), G (E) 
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where Go is the Green operator of the free sp motion. 
Inserting (fT]) with ([2]) into ([J}, and multiplying from the 
left by (s^|t/ _1 (f2 s )e lRflfc , the following set of linear equa- 
tions for the amplitudes A sfJj 
obtained: 
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where H ss > = R s — R s ' . Here there is no direct overlap 
between the two (non-orthogonal) set of basis functions, 



unlike in the secular matrix equation that results from 
the wave function expansion method 14]. Only the off- 



diagonal block (s ^ s') of the matrix for the linear system 
((5|) contains the dependence on the nuclei orientation £! s , 
as U~ 1 (£l s )U(£l s ) is the unitary operator. The matrix 
elements in ([5]) involving the Green operator Go are ex- 
pressed in terms of the Wigner D-functions as follows 

D£ im (n s ) (snljm 1 \G (E)e iR ^' k \s'n , l , j , m2) 

m\rri2 

D'L„A^), (6) 

where the new matrix elements in (j6)) are explicitly given 
in Ref. [13] [see expressions (12)-(16)]. The off-diagonal 
term (s ^ s') of these new matrix elements vanishes 
for large separations [l6j |. turning the two-center prob- 
lem into two independent one-center problems associated 
with the individual deformed nuclei. 

The system of algebraic equations §5§ is equivalent to 
the Lippmann-Schwinger equation (j4|) with separable po- 
tentials ([2|), whose solution is exact. The solvability con- 
dition is that its determinant vanishes, leading to the 
adiabatic energies E that are a parameter in the Green 
operator Go- With the eigenvalues E, the eigenstates \ip) 
are obtained solving the system (0 for the amplitudes 
A SIM and requiring the normalization of the state vectors 
\<p). For well-separated nuclei, the eigenstates are those 
of the individual nuclei, being the projection of the sp 
total angular momentum along their intrinsic symmetry 
axis a good quantum number. Molecular orbitals develop 
at small separations, which may not have good quantum 
numbers (there are no symmetries) as the two set of ba- 
sis states are completely mixed by the potential and the 
operator of finite rotations. For identical potentials (mu- 
tually aligned identical nuclei) , a symmetry of the two- 
center potential {T]) arises, namely its invariance under 
the permutation of the individual orientated potentials 
with respect to the CM point in FigJT] In this case, the 
parity of the molecular sp states is also a good quantum 
number. Symmetric and anti-symmetric linear combina- 
tions of those (asymptotically degenerated) states, with 
opposite parity, result in (atomic) states localized around 
one of the potentials [lfij]. Where this is not the case, the 
solution of the two-center problem at large separations 
directly yields states of the individual nuclei. 

The formalism is now illustrated with calculations of 
molecular sp levels diagram (i.e., energy levels as a func- 
tion of the internuclear distance R) for 12 C + 12 C — > 
24 Mg, where 12 C is an oblately deformed nucleus with 
quadrupole deformation f3 2 = —0.5 [23]. Calculations 
are carried out for the collisions shown in Figj2] The pa- 
rameters of the asymptotic WS potentials (R = and 
R — » oo) including the spin-orbit term are the global 
parameters by Soloviev [2J], where the depth of the po- 
tentials has been adjusted to reproduce the experimental 
sp separation energies [25[ . To describe fusion, all the po- 
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FIG. 2: Different configurations for the reaction 12 C + 
axial symmetric (a), and non-axial symmetric (b)-(d). 
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FIG. 3: (Color online) Neutron molecular sp levels as a func- 
tion of the internuclear distance R, corresponding to the con- 
figurations shown in FigO Different curves in panel (a) are 
associated with different magnetic quantum numbers. The 
points denote the Fermi level of the spherical compound nu- 
cleus and the 12 C nuclei. 



tential parameters (including those of the Coulomb po- 
tential for protons) have to be interpolated between their 
values for the separated nuclei and the spherical com- 
pound nucleus [T7(. The parameters can be correlated 
by conserving the volume enclosed by a certain equipo- 
tential surface (the Fermi level of the spherical fused sys- 
tem) of the two-center potential (fl]), for all separations 



R between the nuclei with orientations f2 s 14. .Irj. The 
nuclear shape is considered to be the same for neutrons 
and protons, with the neck size naturally determined by 
the superposition of the smooth tail of the two WS po- 
tentials. However, a variable neck size could be included 
in the method through an additional intermediate po- 
tential, resulting in a three-center problem. This may be 
very useful in fission studies and cluster physics. 
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FIG. 4: (Color online) The same as Fig(3] but for protons. 



Figures [3] and 0] show the molecular sp levels diagram 
for neutrons and protons, respectively, for different ori- 
entations of the two oblately deformed 12 C, as presented 
in FigO Only for the aligned orientation of the defor- 
mation axis with the internuclear axis [axial symmetric 
configuration in panel (a)], the projection of the nucleon 
total angular momentum along the internuclear axis is 
a good quantum number at all separations, whose val- 
ues are represented by different curves. Here as in panel 
(b) , the parity of the molecular orbitals (not indicated) 
is also a good quantum number, whose values are eas- 
ily deduced from the orbital angular momentum I of the 
sp levels of the spherical compound nucleus (R = 0), 
i.e., (—I) 1 - The full circles denote the Fermi level of the 
12 C nuclei and the spherical compound nucleus. As ex- 
pected, the asymptotic shell structure does not depend 
on the mutual alignment of the 12 C nuclei. 

These molecular spectra show significant features: 

(i) The asymptotic shell structure of 12 C is much less 
distorted for non-axial symmetric configurations 
[panels (b)-(d)] than for the axial symmetric one 
[panel (a)], the former showing more preservation 
of the indentity of the overlapping nuclei. 

(ii) For non-axial symmetric configurations, the neu- 
tron levels show a minimum at separations between 
4-6 fm, which may result in a "molecular pocket" 
in the collective potential energy surface [8[. This 
depends weakly on the mutual alignment. 

(iii) Many avoided crossings appear between 2-6 fm, in 
which the sp wave-function abruptely changes its 
nodal structure. It may lead to strong peaks in 
the radial collective mass parameter [27| , which can 
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hinder the fusion of the nuclei. The critical radius 
for fusion [26|] (where the shell structure of 24 Mg 
starts to build up) is quite small for non-axial sym- 
metric configurations (< 2 fm), which also favours 
the formation of a nuclear molecule (the weakly 
overlapping nuclei remain longer around 4-6 fm). 

The preservation of the identity of the nuclei along 
with their being trapped around the contact distance are 
crucial aspects for the formation of a nuclear molecule 
Q. These favourable features are shown by non-axial 
symmetric configurations of 12 C + 12 C. Since their shell 
structures are quite similar, the 12 C orientation is clearly 
an essential variable in the reaction processes. This de- 
gree of freedom activates 28|, 29} dynamical modes (but- 
terfly, anti-butterfly, belly-dancer etc) at contact, making 
the nuclei "dance" there for quite some time. It should re- 
sult in narrow resonances in the reaction cross sections, as 
shown by measurements [3(]} . The molecular sp spectra 
are useful to microscopically obtain collective potentials 
and mass surfaces for a molecular reaction dynamical cal- 
culation [29|1 • It is worth mentioning that other studies 




have also argued the importance of the 
non-axial symmetric configurations for effects of molecu- 
lar resonances on reaction processes for 12 C + 12 C. 

In summary, a general new technique to solve the two- 
center problem with arbitrarily-orientated deformed re- 
alistic potentials has been demonstrated. Among other 
applications such as in cluster physics [35} , this should be 
very useful for describing, within the molecular picture, 
low energy nuclear reaction processes involving deformed 
nuclei, such as (i) formation of heavy and superheavy el- 
ements [111 i (ii)effects of breakup of weakly-bound nu- 
clei on fusion [37]| . and fusion reactions of great astro- 
physical interest [27|. Molecular sp spectra clearly show 
that non-axial symmetric configurations are crucial for 
the formation of a nuclear molecule in the reaction 12 C 
+ 12 C. Reaction dynamical calculations for a quantita- 
tive understanding of molecular resonance structures in 
its astrophysical S-factor [30( are in progress. 

The author thanks Prof. W. Scheid for discussions. 
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[1] J.E. Howard et al., Phys. Rev. A 52, 4471 (1995). 
[2] J.S. Slater, Electronic Structure of Molecules, Mc- 
GrawHill, New York, 1963. 



[3 

[4 

[5. 
[6 

[-: 

[8 

[9 
[10 

[11 
[12. 

[13 
[14 

[15 



[16 
[17 

[is: 

[19 

[20 
[21 

[22' 
[23 
[24 

[25' 
[26 
[27 

[28 

[29 

[30 
[31 
[32 
[33 
[34 

[35 

[36 

[37 



B. Miiller, J. Rafelski and W. Greiner, Phys. Lett. B 47, 
5 (1973). 

S. Chandrasekbar, Proc. R. Soc. London Ser. A 421, 227 
(1989). 

E. Teller, Z. Phys. 61, 458 (1930). 
J. Eichler, Phys. Rev. Lett. 75, 3653 (1995). 
Da-Heng He et al., Phys. Rev. D 70, 094004 (2004). 
W. Greiner, J.Y. Park and W. Scheid, Nuclear Molecules, 
World Scientific, Singapore, 1994. 

J.A. Maruhn and W. Greiner, Z. Phys. 251, 431 (1972). 
M. Mirea, Nucl. Phys. A 780, 13 (2006). 
RA. Gherghescu, Phys. Rev. C 67, 014309 (2003). 
Yu.Ts. Oganessian, Nature 413, 122 (2001). 

D. J. Hinde et al., Phys. Rev. Lett. 100, 202701 (2008). 
G. Nuhn, W. Scheid and J.Y. Park, Phys. Rev. C 35, 
2146 (1987). 

J. Revai, JINR, E4-9429, Dubna, 1975; B. Gyarmati, 
A.T. Kruppa and J. Revai, Nucl. Phys. A 326, 119 
(1979). 

FA. Gareev et al., Nucl. Phys. A 286, 512 (1977). 

A. Diaz- Torres and W. Scheid, Nucl. Phys. A 757, 373 

(2005). 

C. E. Rolfs and W.S. Rodney, Cauldrons in the Cosmos, 
The University of Chicago Press, Chicago, 1988. 

A. Faessler and R.K. Sheline, Phys. Rev 148, 1003 
(1966). 

T. Tamura, Rev. Mod. Phys. 37, 679 (1965). 

E. Merzbacher, Quantum Mechanics (Second Edition), 
John Wiley and Sons, Inc., New York, 1970, p. 371 

A. Diaz- Torres, Phys. Lett. B 594, 69 (2004). 
V.M. Lebedev et al., Phys. Atom. Nuclei 62, 1455 (1999). 
V.G. Soloviev, Theory of Complex Nuclei, Pergamon 
Press, Oxford, 1976, p. 21 

G. Audi and A.H. Wapstra, Nucl. Phys. A 565, 1 (1993). 

D. Glas and U. Mosel, Phys. Lett. B 49, 301 (1974). 

A. Diaz-Torres, L.R. Gasques and M. Wiescher, Phys. 
Lett. B 652, 255 (2007). 

P.O. Hess, W. Greiner and W.T. Pinkston, Phys. Rev. 
Lett. 53, 1535 (1984). 

E. Uegaki and Y. Abe, Phys. Lett. B 231, 28 (1989); 
Phys. Lett. B 340, 143 (1994). 

E.F. Aguilera et al., Phys. Rev. C 73, 064601 (2006). 

H. Chandra and U. Mosel, Nucl. Phys. A 298, 151 (1978). 
O. Tanimura and T. Tazawa, Phys. Lett. B 78, 1 (1978). 
M. Niklas et al., Z. Phys. A 323, 27 (1986). 

P.O. Hess, J. Schmidt and W. Scheid, Ann. Phys. (NY.) 
240, 22 (1995). 

W. von Oertzen, M. Freer and Y. Kanada-En'yo, Phys. 
Rep. 432, 43 (2006). 

A. Diaz- Torres, Phys. Rev. C 69, R021603 (2004); Phys. 
Rev. C 74, 064601 (2006). 

A. Diaz- Torres et al., Phys. Lett. B 533, 265 (2002); 
Nucl. Phys. A 703, 83 (2002). 



